Accepted by ApJ, 23 January 2013 

Preprint typeset using LAT^X style cmulateapj v. 5/2/11 



GBT DETECTION OF POLARIZATION-DEPENDENT HI ABSORPTION AND HI OUTFLOWS IN LOCAL 

ULIRGS AND QUASARS 

Stacy H. Teng 1- 7 , Sylvain Veilleux 2 ' 3 ' 4 ' ''■ 7 , & Andrew J. Baker 6 

(Received; Accepted) 
Accepted by ApJ, 23 January 2013 

ABSTRACT 

We present the results of a 21-cm HI survey of 27 local massive gas-rich late-stage mergers and merger 
remnants with the Robert C. Byrd Green Bank Telescope (GBT). These remnants were selected from 
the Quasar/ULIRG Evolution Study (QUEST) sample of ultraluminous infrared galaxies (ULIRGs; 
L8-1000 urn > 10 12 Lq) and quasars; our targets are all bolometrically dominated by active galactic 
nuclei (AGN) and sample the later phases of the proposed ULIRG-to-quasar evolutionary sequence. 
We find the prevalence of HI absorption (emission) to be 100% (29%) in ULIRGs with HI detections, 
100% (88%) in FIR-strong quasars, and 63% (100%) in FIR-weak quasars. The absorption features 
are associated with powerful neutral outflows that change from being mainly driven by star formation 
in ULIRGs to being driven by the AGN in the quasars. These outflows have velocities that exceed 
1500 km s _1 in some cases. Unexpectedly, we find polarization-dependent HI absorption in 57% of our 
spectra (88% and 63% of the FIR-strong and FIR-weak quasars, respectively). We attribute this result 
to absorption of polarized continuum emission from these sources by foreground HI clouds. About 
60% of the quasars displaying polarized spectra are radio-loud, far higher than the ~10% observed 
in the general AGN population. This discrepancy suggests that radio jets play an important role in 
shaping the environments in these galaxies. These systems may represent a transition phase in the 
evolution of gas-rich mergers into "mature" radio galaxies. 

Subject headings: galaxies: active — galaxies: evolution — quasars: absorption lines — ISM: jets and 
outflows — radio lines: galaxies 



1. INTRODUCTION 

Recent theoretical work suggests that feedback from 
AGN plays an important role in establishing the present- 
day appearance of massive galaxies. Major galaxy merg- 
ers at high redshifts may trigger powerful starbursts, 
lead to the formation of elliptical galaxies, and ac- 
count for the growth of supermassiye black holes (e.g., 
Sanders et "aJ]ll988Hdi Matteo et aIll2005l:[Springel et alJ 
20051: iHopkins et alJ 120051 120081 ). In this merger-driven 
evolutionary scenario, two gas-rich progenitors collide to 
form a completely obscured ULIRG where substantial in- 
flows of gas trigger starbursts and feed the central black 
hole(s). As the system evolves, the obscuring gas and 
dust begin to disperse, bringing about the formation of a 
dusty IR-excess quasar and eventually a "naked" optical 
quasar. Powerful winds, driven by either the starburst 
or the central black hole, seem to be required to quench 
star formation and to stop the growth of the black hole 
in order to explain the observed Mbh — & relation (e.g., 
IMurravet afllSoOfil L 
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There is growing support for the existence of this type 
of "negative feedback" in luminous late-stage mergers. 
Most galaxies with star formation rate (SFR) densi- 
ties >0.1 Mp) yr" 1 kpc~ 2 show sig ns of outflow (e.g., 
iHeckmanl l2002t iVeilleux et all [2005) . These starburst- 
driven winds are often seen in ULIRGs, with mass out- 
flow rates of ~0.1-5 times the SFR and vel ocities ~100- 
400 km s" 1 (IRunke et al. 1 [20021 l27)05aM l IMartinl l2(M 
2006). In late-stage mergers and pure Seyfert galaxies, 
however, these winds are sometimes driven by the AGN, 
reaching outflow velocit ies >1000 km s" 1 (Rup ke~et al.1 
[2005cl IKrug et al.ll2f)Toh . Similar AGN-driven outflows 
are also found in recently (re)started radio gal axies (e.g., 
lOosterloo et~aT1 120001 iMorganti et "all [2005b). Several 
studies have also shown that these outflows are of- 
ten multi-phased, with powerful molecular, neutral, and 
ionized outflo ws having all been detected in the same 
object (e.g.. IMorganti et all l2007t iFischer et al.1 120101: 
Feruglio et al . 2010 HRupke fc Veilleuxll2011l : ISturm et al.1 

201 in. 

In the present paper, we focus on the HI properties of 
the QUEST ULIRGs and quasars, searching for the pres- 
ence and the driving source (AGN or starburst) of neu- 
tral outflows in merging systems as a function of merger 
stage. This paper is organized as follows. In Section [2j 
we discuss how our HI sample is selected. In Section [3l 
we describe our observational setup, data reduction, and 
spectral analysis. We then examine the HI profiles of 
our targets in Section [4] and the detection of polarized 
absorption (including a more sophisticated spectral de- 
composition that boosts our number of absorption de- 
tections) in Section [5j Finally, in Section [6j we discuss 
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the possible origins of the HI outflows. We conclude in 
Section [7] with a summary of our findings. Throughout 
this paper, we adopt H n = 71 km s" 1 M pc -1 , £Im — 
0.27, and fl A = 0.73 (jHinshaw et al.ll2009h . 

2. SAMPLE SELECTION 

In order to assess the physical conditions and the dy- 
namical evolution of the neutral gas in late-stage merger 
systems, we selected targets from the QUEST sample 
of ULIRGs and quasars that are bolometrically domi- 
nated by AGN. The details of the QUEST sample se- 
lection are outlined i n previous Q UES T papers (e.g., 
Schweitzer et all 120061: fNetzer et alj|2007t iVeilleux et al.l 
2009aD. Briefly, the original QUEST sample of ULIRGs 
begins with the 1-Jy sample, a flux-limited sample of 
118 ULIRGs selected from a redshift survey of the IRAS 
faint source catalog (|Kim fc Sanders! Il998| ). To these, 
five well-studied local ULIRGs from the Re vised Bright 
Galaxy Survey (RBGS; ISanders et afll200"l were added 
given the availabi lity of archival data from multiple 
observatories (see IVeilleux et all I2009at IVeilleux! 120121 
for details). The original quasar sample contains 25 
quasars, including 24 Palo mar-Green (PG) quasars from 
the Bright Quasar Sample (Sch midt fc G reen 1983]) and 
B2 22 01+31 A which satisfies the B magnitude crite- 
rion of Schmid t fc Green] (fl98 3). An additional nine PG 
quasars were added to the sample from archival data. 
The selected quasars fall at the low end of the quasar 
B-band luminosity function. 

The ULIRGs and quasars in the QUEST sample are 
well-matched in redshift (z < 0.3) and in near-IR imag- 
ing show tidal features th at indicate merger origins 
(jVeilleux et al.ll2?M I2009bh . The entire QUEST sam- 
ple spans the merger evolution sequence from binary 
starburst-dominated ULIRGs, to completely coalesced 
AGN-dominated ULIRGs, to dust-enshrouded quasars, 
to optical quasars whose dust has been cleared. 

Our HI sample is a subset of the QUEST targets. To 
ensure quality spectra, we require z < 0.12; our 27 tar- 
gets evenly sample coalesced ULIRGf0 (9 targets), far- 
infrared(FIR)-strong quasars (/60/xm//i5/xm > 1; 10 tar- 
gets) , and FIR- weak quasars (8 targets) . The two classes 
of quasars have similar spectral energy distributions rep- 
resentative of AGN, with the FIR-strong quasars having 
an additional contribut ion to their FIR e mission from 
intense star formation (jNetzer et al.l l2007h . The mor- 
phologies of the quasar hosts are in agreement with the 
picture that FIR-weak quasars are at a later evolution- 
ary stage than FIR-strong quasars, following th e (near) 
cessation of star formation ()Veilleux et al.|[2009bf) . Thus, 
our targets sample the later stages of merger evolution 
from coalescence to the aftermath of significant star for- 
mation. Table [T] summarizes the properties of our HI 
sample. 

3. OBSERVATIONS AND DATA 

3.1. Observational Setup 

The targets were observed with the GBT Spectrometer 
in the L band between 2011 August and 2012 March dur- 
ing the ^160 hours of time allocated for GBT programs 
11B075 and 12A075 (PI: Teng). The data were collected 

8 This category includes Mrk 273, the only source in our sample 
to have binary nuclei (see Appendix). 



in 10 minute on-off source pairs (5 minutes per posi- 
tion) with 12.5 MHz bandwidth, two spectral windows 
centered at the same frequency (1420.4058 MHz), two 
linear polarization modes (XX and YY), and nine sam- 
pling levels. To minimize the effects of radio frequency 
interference (RFI) on the data, we used the notch filter 
for sources with heliocentric recessional velocities below 
17000 km s _1 . The total on-source time for each galaxy 
was determined during the observing runs based on the 
strength of the spectral features relative to the noise or 
up to a maximum of three hours. Most sources were ob- 
served over multiple observing sessions. The circumpolar 
starburst galaxy M82 was observed at the beginning of 
each observing session as a check for instrumental varia- 
tions. 

We also measured the continuum flux densities of our 
sample at 1.42 GHz using the AutoPeak procedure in 
order to check for time variability. However, due to 
the plethora of continuum sources within the large GBT 
beam, the 5-cr confusion limit is ^100 mJy at this fre- 
quency. Therefore, we were only able to obtain confi- 
dent continuum measurements for five targets. These 
values are in agreement with older values obtained with 
the Very Large Array (VLA) for their nuclear core flux 
densities (|Condon et al.l [20021) . despite the much larger 
beam of the GBT. For the remainder of sample, we use 
published values of the core con tinuum flux densities 
as measured with the VLA (e.g . . | Barvainis et al.1 119961 : 
ICondon et alJH998t IRafter et al.ll2009fi ~ 

3.2. Data Reduction 

The spectra wer e reduced using GBTIDL 
(jMarganian et all 12006') . Individual five-second records 
showing large broadband harmonic RFI were flagged, 
and persistent RFI spikes that appear in the data were 
replaced by interpolation across the affected channels 
for both spectral windows and polarizations. For each 
polarization, the data from the two spectral windows 
were averaged together, and a polynomial of order <5 
was fitted over a range of ^1000 km s _1 on each side of a 
galaxy's systemic redshift in the Hanning-smoothed and 
then decimatec0 spectrum to subtract the baseline. The 
same baseline-fitting regions and order of polynomial fit 
were used for both polarizations. 

Flux calibration was derived from simultaneous ob- 
servations of an internal noise diode. The baseline- 
subtracted and flux-calibrated scans for each polariza- 
tion were accumulated and averaged for each galaxy. A 
fourth order boxcar smoothing was then applied to the 
averaged data with decimation, resulting in a resolution 
of ^6 km s _1 per channel. The two polarizations were 
then averaged together to produce the full intensity spec- 
tra shown in Figure [T] We ensured that our flux calibra- 
tion is correct by comparing the full intensit y profile of 
Mrk 231 with that measured by the VLA (jCarilli et al.l 
1998); both profiles are identical. Discussions of individ- 
ual objects in the Appendix include comparisons with 
previous data, whose generally good agreement with our 
data in terms of velocity profile and flux validate our 
data reduction process. Table [1] also provides a list of 
some observational parameters for our galaxies. 

9 The total number of channels in the spectrum is decreased, or 
decimated, by a factor equal to the order of smoothing performed. 
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Nearly all of the galaxies were detected. Two galax- 
ies, PG 0804+761 and F15462-0450, had strong RFI 
spikes overlapping with the expected HI frequencies. At 
the limit of three hours on source, PG 1126-041 and 
F21219-1757 were not detected, and F04103-2838 and 
F15130-1958 were marginally detected. The detections 
of PG 1426+015 and PG 1617+175 are uncertain be- 
cause their putative emission features have widths sim- 
ilar to those of broadness of the baseline ripples before 
baseline subtraction. 

3.3. Spectral Analysis 

The individual HI profiles in the spectra of our sam- 
ple galaxy are complex combinations of absorption and 
emission. For simplicity, we have measured the emis- 
sion and absorption features separately, but we note that 
this approach may slightly alter the final measured val- 
ues of these quantities. More complex decomposition of 
the emission and absorption is only possible for a few 
sources because of the degeneracy between emission and 
absorption, as discussed in Section 1531 

For galaxies with significant emission (peak-to-rms > 
3a), we measure the central velocity (Vq), width (W20, 
W50), and total velocity-integrated flux (/hi) using the 
GBTIDL task gmeasure. In these instances, since there 
are no obvious companions enclosed within the GBT 
beam, the total flux is integrated between velocities at 
which the flux density reaches zero. The central velocity 
is the heliocentric velocity measured at the midpoint of 
the 20% flux level. The width of emission is measured 
at the 20% and 50% levels (W20 and W50, respectively). 
The measurement errors in these quantities a re der ived 
using the expressions prese nted inlWei et al.l (2010) fol- 
lowi ng the method o logies of lSchneider et al.l (|1986l . ll990| ) 
and iFouque et al.l (|1990f ) . The HI c olumn density is 
calcul ated using expression 12.17 from lRohlfs fc Wilson! 
(2000) for optically thin radiation. Assuming the GBT 
antenna temperature to flux density gain conversion of 
-2 K Jy" 1 : 

Am = 3.64 x 10 18 /hi cm- 2 , (1) 

where /hi is the velocity integrated HI flux in units of 
Jy km s" 1 . The neutral ga s mass is calculated following 
IHavnes fc Giovanellil (Il98l : 

M m = 2.35 x 10 5 /hi(-^) 2 M o , (2) 
tin 

where V cor is the recessional velocity corrected for the 
influence of the Virgo Cluster, the Great Attractor, and 
the Shapley Supercluster as calculated by the NASA Ex- 
tragalactic Database (NED). All of these measured and 
derived quantities are listed in Table [2l 

Only a few absorption features are significantly de- 
tected (by the criterion that peak-to-rms > 3a). How- 
ever, many of the absorption features are very broad but 
shallow, so the peak-to-rms ratio is not necessarily the 
best determinant of a detection. We also consider a fea- 
ture to be detected if its full width at half maximum 
(FWHM) spans at least 10 velocity channels and it has a 
peak-to-rms > la. These absorption features are quan- 
tified using the GBTIDL task fitgauss, with one Gaus- 
sian component at a time added to a spectrum until the 
residuals are statistically featureless. With the excep- 



tion of PG 1119+120, where a narrow absorption feature 
is present within the emission profile, only the absorp- 
tion features are modeled. For PG 1119+120, the emis- 
sion was first modeled with Gaussians and subtracted 
from the data to define a baseline. The absorption was 
then modeled from the emission-subtracted spectrum. 
The Gaussian models are plotted (as components and 
as sums) in Figure [TJ The GBTIDL task measures the 
central velocity (V Q ), FWHM, the "peak" flux density, 
and their measurement errors for each Gaussian compo- 
nent. The full velocity-integrated absorbed flux (/ a bs), 
optical depth (t), and column density (Nm) are then 
calculated from the best-fit models. We assume the HI 
absorbing screen covers the nuclear continuum source, in 
which case the optical depth profile is: 

r(t,) = -Ml-^L], (3) 

J1A GHz 

where Sin(v) is the HI line profile as a function of 
velocity and 6*1.4 GHz is t he continuum flux density 
( Whitcoak & Gardner 1972). The column density is then 
calculated as: 

^ = 1.823 x 10 18 f t(v) dv cm" 2 , (4) 

1 s J 

where T s is the spin temperature. In this paper, we will 
assume a spin temperature of 1000 K (see Section IB - !]) . 
For absorption features in PG 1119+120, PG 1244+026, 
and PG 1440+356, the absorption troughs are deeper 
than their nuclear continua, suggesting that sources of 
continuum emission extend beyond the nucleus. Table [3] 
summarizes the measured and derived values for all of 
the detected absorption features. 

As a reality check, we inspected the data to confirm 
that the choice of polynomial order for the baseline fitting 
does not affect the detection of absorption or emission. 
Specifically, we examined whether polynomial orders of 3 
or 5 can introduce ripples into the data that masquerade 
as absorption or emission features. After repeating the 
baseline subtraction for our data using a range of polyno- 
mial orders (1-5), we confirm that the same spectral fea- 
tures are still detected, even for marginal detections such 
as that in F04103-2838. The measured widths and veloc- 
ities are in agreement with those measured with the best- 
fit baseline to within the errors. Therefore, we are confi- 
dent that even the weak detections noted as "marginal" 
in Tables [5] and [3] are real. 

4. HI PROFILES OF LATE-STAGE MERGERS 

The individual HI profiles of our targets are discussed 
in the Appendix. In general, we see HI trends along the 
merger sequence but no correlation with infrared lumi- 
nosity. 

All of the detected ULIRGs show HI absorption, typi- 
cally with multiple velocity components, many of which 
are narrow (FWHM < 100 km s" 1 ) and some of which 
are redshifted. Only 2/7 (29%) of the detected ULIRGs 
have HI emission. The different velocity absorption com- 
ponents have different origins and are not necessarily 
caused by outflows or infall; kinematics of circumnuclear 
gas disks may mimic the properties of outflows or in- 
fall in single dish data. In a couple of cases, such as 
UGC 05101 and Mrk 273, the absorption profiles are 
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nearly symmetrical about zero velocity. For the Type 2 
AGN UGC 05101, the profile is likely from absorption 
due to a slightly inclined rotating HI disk superposed 
on an obscured nucleus. This interpretation is consis- 
tent with CO maps of UGC 05101 s howing smooth, reg- 
ular rotation of th e molecular gas (jGenzel et al.1 119981 : 
I Wilson et aTll2008[ ). On the other h and, the CO velocity 
field of Mrk 273 is very distorted (|Downes fc Solomon! 
Il998t iWilson et ai1l2008[ ). so it is unlikely that the gas 
is in well-ordered disk-like rotation. The double absorp- 
tion in HI seen in Mrk 273 probably corresponds to the 
binary obscured nuclei in this object. Thus, CO maps 
can help eliminate gas disks as the origin of the absorp- 
tion features. However, only a few galaxies in our sample 
have readily available CO maps so direct comparisons are 
difficult. Comparison of single dish CO and HI profiles 
for the QUEST targets is on-going, but it is beyond the 
scope of the present paper. The more complex profiles, 
such as those seen in F15250+3608, F05189-2524, and 
F15130-1938, may be interpreted as multiple infalling or 
outflowing components. They may also result from the 
superposition of an emission profile and a broad absorp- 
tion profile. The latter case is discussed in more detail 
in Section l5~2l 

In contrast, the HI profiles of the FIR-strong quasars 
are generally detected not only in absorption (88% or 7 of 
the 8 detected sources) but also in emission (75%). The 
absorption features are on average —60 km s _1 broader 
than those seen in ULIRGs or approximately three times 
greater than the typical FWHM measurement error in 
the ULIRGs. Both redshifted and blueshifted absorp- 
tion features are seen, although the systemic velocities 
of these systems are uncertain. Optical emission lines 
may underestimate recessional velocity since some frac- 
tion of the line emission may arise from outflowing ma- 
terial. In the case of PG 0050+124 (I Zw 1), this effect 
causes a shift of —700 km s _1 (see Appendix). A simi- 
lar situation could also hold for PG 1440+356. In three 
cases, PG 1119+120, PG 1244+026, and PG 1440+356, 
the minima of the absorption troughs exceed the nu- 
clear continuum flux densities. This discrepancy implies 
that off-nuclear emission, not accounted for in the nu- 
clear fluxes of Table [TJ also contributes to the contin- 
uum emission in these four sources. Given resolution 
limit of the VLA core measurements of PG 1119+120 
and PG 1440+356, the off-nuclear emission must be at 
least -20-30 kpc from the core; for PG 1244+026, the 
separation from the core must be at least 3 kpc. While 
we cannot formally rule out the possibility that these ab- 
sorption features are absorption of HI against continuum 
from background quasars, this possibility is statistically 
unlikely. Another possibility is that the flux densities of 
the cores in these galaxies are time variable as exempli- 
fied by the variability of the core in PG 1351+640 (see 
Appendix) . 

The HI profiles of the FIR-weak quasars are different 
from those of the ULIRGs and FIR-strong quasars. HI is 
seen mostly in emission (88% or 7 of 8 detected sources), 
and only 25% (2/8) show broad absorption. The emis- 
sion profiles are very broad (IU50 >500 km s _1 ) and 
bumpy, perhaps beginning to resemble those of regularly 
rotating disks that typically display Gaussian or double- 
horned profiles; the measured widths and velocities are 



consistent with previous measurements of quasars (e.g., 
Bothun et aLlll984: iCondon et al.l^985t iHutchi ngs et al.1 
19871 IHq et al.ll2008t iKonig et al.l 120091 ) . These proper- 
ties will be revisited in Section 1531 

Considering now our entire HI sample, the distribu- 
tion of equivalent widths (EW; defined as the integrated 
flux density for each component divided by the contin- 
uum flux density of the core) is plotted in Figure [2] For 
this definition, the equivalent widths of the absorption 
components are negative. It is clear that the strengths 
of the emission and absorption features progress along 
the merger sequence, with ULIRGs having the weakest 
features, the FIR-strong quasars having the strongest 
absorption features, and the FIR-weak quasars having 
the strongest emission features. In Figure |3l we com- 
pare the FWHMs of emission and absorption with the 
star formation rates of their host galaxies. For emis- 
sion, the FIR-weak sources have the broadest widths, 
the ULIRGs have the narrowest widths, and the FIR- 
strong quasars bridge the gap between those two classes. 
In the ULIRGs, there appears to be a weak correlation 
between the width of the absorption components and the 
SFR, reminiscent of a similar relationship seen in Na I 
(|Rupke et al.ll2005bf ). However, given that multiple ab- 
sorption components are associated with a single system 
(and thus a single SFR value), the validity of the connec- 
tion between the SFR and the FWHM of the absorption 
components in ULIRGs is uncertain. On the other hand, 
absorption FWHM is not at all related to SFR in the 
quasars, further suggesting that the AGN play a role in 
these systems. 

5. POLARIZATION-DEPENDENT HI ABSORPTION 

An unexpected outcome of our GBT program is that 
many of the absorption features show a dependence on 
polarization (Figure |4j. Indeed, 13 of our 23 GBT- 
detected targets show some profile difference in their two 
polarizations (1 ULIRG, 7 IR-excess quasars, and 5 opti- 
cal quasars). For GBT observations, it is expected that 
there is a <6% difference between polarizations in terms 
of normalization, but the observed spectral shapes should 
be the same in both polarizations (R. Maddalena, private 
comm.). This expectation is confirmed by the fact that 
our M82 data remain in agreement after baseline sub- 
traction. Polarization variation in our QUEST sample 
is not RFI-related, as we have removed time intervals 
where the data are affected by broadband RFI. Since it 
is improbable that the HI screen itself is polarized, the 
underlying continuum emission must be polarized. 

5.1. Time Variability 

The most dramatic examples of polarized spectra are 
seen in PG 1211+143 and PG 1440+356. Not only are 
these spectra polarized, but they are also time variable 
on timescales of —30 minutes (Figures [5] and [5]) ■ In fact, 
of our radio-detected sample galaxies, time variability 
is only seen in sources that have significant absorption 
and spectra that are significantly polarized. A possi- 
ble concern is that since we observed these sources over 
several sessions, variability could be an instrumental ef- 
fect. However, comparing the different sessions' data for 
M82, we find that time variability and polarization are 
not seen beyond what is expected of typical GBT obser- 
vations. Moreover, session-to-session instrumental vari- 
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ability would not explain the 30-minute variability seen 
within the same observing session. Time varying polar- 
izations are not seen only in the targets with the brightest 
continua, so the variations are not likely to be due to in- 
strumental systematics from poor cance llation of base line 
ripples in the position-switching scans (|O'Neilll2002K 

A different sort of instrumental effect may be leakage 
of linear polarization from the linearly polarized back- 
ground continuum into the total intensity of the absorp- 
tion profiles. We cannot yet exclude this scenario since 
these observations were not designed for polarimetry and 
our current data lack the cross-correlated products (XY 
and YX) that would be needed to confirm these time 
variations. 

If these time variations are real, they could be caused 
by one of several processes. First, they could be due to 
variations in the compact continuum sources behind the 
HI screen. As the X-ray light curves of PG 1211+143 
and PG 1440+356 show (Figure 0), the X-ray emitting 
central region can be highly variable on a timescale of 
about 30 minutes. Although not contemporaneous with 
our GBT data, these X-ray data suggest that the AGN 
continuum can vary on a 30-minute time scale if the radio 
continuum is indeed related to the X-ray producing re- 
gions. However, if the polarization-dependent variations 
are simply a result of the changing compact source con- 
tinuum, then we would expect the continuum to vary in 
block and would observe a change in the equivalent width 
of the HI absorption, not the frequency-dependent vari- 
ations we see in these targets. Second, variability could 
be a result of changes in the structure of the HI screen. 
However, the timescale for the destruction of small HI 
clouds in the interstellar medium is on the order of a 
Myr, depending on the exact propertie s of the interstel- 
lar medium causing the absor ption (e.g. JMcKee fc Cowid 
H977tlNagashima et aD2006[ ): it is therefore unlikely that 
cloud destruction causes the observed short-term varia- 
tions. Similarly, it is unlikely that the movement of HI 
clouds into and out of our line of sight to the continuum 
source can cause the very broad variations observed in 
these objects. 

The time variability of the polarization-dependent HI 
absorption is easier to explain if the polarized contin- 
uum emission comes from both the nuclear source and an 
extended source (e.g., radio jets). In this picture (Fig- 
ured]), the HI screen is also extended, covering both the 
nuclear source and extended jets with kinematically dis- 
tinct clouds. As the central source fluctuates, the ratio 
of central-to-jet emission also changes, resulting in the 
observed wavelength-dependent polarization variations. 
In this scenario, there is no need for the structure of the 
HI screen to change on a fast timescale. 

5.2. Decomposing the Spectrum 

Regardless of whether the time variability is real, it is 
certainly the case that the spectra shown in Figure U] are 
polarized. Again, we checked for baseline-subtraction- 
related effects on the detection of polarized absorption 
and found that the choice of the polynomial order of 
the baseline fit matters little to the detection of these 
polarization-dependent absorption features. The most 
obvious example of polarization not due to observational 
systematics is seen in the spectra of PG 0050+124 (I Zw 
1). This HI profile is not polarized at relative velocities 



of 500-1000 km s _1 but becomes polarization-dependent 
between relative velocities -500 and +500 km s _1 (Fig- 
ure [4]) . This behavior argues against attribution of ap- 
parent polarization dependencies to poor baseline sub- 
traction or session-to-session variations in our GBT ob- 
servations. At least in this one case, the polarization de- 
pendence must be due to something physically happen- 
ing at relative velocities between -500 and +500 km s -1 . 
When we compare the HI profiles of each galaxy in the 
two separate polarizations, it appears that one polariza- 
tion typically has higher flux than the other, with the 
exception of PG 1440+356, in which there is a veloc- 
ity shift instead. The effect is not systematic in the 
sense that one polarization is always higher than the 
other; each is nearly equally likely to be the polarization 
with the higher intensity. In some cases, particularly in 
PG 1351+640 and PG 1411+442, the lower-intensity pro- 
file appears to be the same as the higher-intensity profile 
with the addition of a broad absorption component. 

To follow up on this intriguing hint, we produced the 
difference spectra (the spectrum in the polarization with 
the more negative HI profile minus the spectrum in the 
polarization with the more positive HI profile) shown 
in Figure |9] The absorption features in these spectra 
are then modeled the same way as before, using the 
GBTIDL task fitgauss. When these absorption mod- 
els are included, many of the full intensity spectra that 
need multiple components to explain the blue- and red- 
shifted components can be explained simply by HI emis- 
sion that is combined with a single broad absorption com- 
ponent. This pattern is most obvious in PG 1351+640 
and PG 1411+442. The results of this type of complex 
decomposition of the spectra are provided in Tables|4]and 
[5l Similar decompositions can explain the shapes of the 
full intensity spectra of the ULIRGs F05189-2524 and 
F15130-1958, and the FIR-strong quasar PG 1244+026. 
However, since those spectra are not polarized, we can- 
not create difference spectra to model the broad absorp- 
tion component. The fitgauss task statistically prefers 
the multiple-component models to the single- absorption 
model for these spectra. Overall, with the exception of 
PG 1501+106, the decompositions result in stronger fea- 
tures in both absorption and emission. 

Given the results of this exercise, we find the instances 
of HI absorption (emission) to be 100% (29%) in detected 
ULIRGs, 100% (88%) in FIR-strong quasars, and 63% 
(100%) in FIR-weak quasars. The more complex de- 
composition of the spectra does not drastically change 
the distribution of EWs (Figure ITUj) compared to Fig- 
ure O Figure [TT] provides a comparison of the FWHMs 
of the emission and absorption features from the complex 
modeling with the star formation rates of these galaxies. 
As in Figure El the FIR-weak quasars have the broad- 
est emission features and the ULIRGs have the narrow- 
est emission features. In terms of absorption, there is 
no obvious relationship between the absorption feature 
widths and the star formation rates. As before, the FIR- 
weak quasars have the broadest absorption features, with 
the FIR-strong quasars spanning the range of widths be- 
tween FIR-weak quasars and ULIRGs. 

6. THE PROPERTIES OF HI OUTFLOWS 

As discussed in Section [4] the absorption components 
are not necessarily caused by infall or outflow. It is pos- 



G 



Teng, Veilleux, & Baker 



sible, as seen in UGC 05101 (edge on) and Mrk 231 (face 
on) , that these features can be associated with the kine- 
matics of gas disks. A circumnuclear gas disk seen nearly 
edge on with a non-isotropic gas mass (or background 
continuum) distribution can mimic characteristics of in- 
fall and outflow. Unfortunately, our single dish data can- 
not unambiguously exclude gas disks as the source of the 
features we observe, but it is unlikely that all of our tar- 
gets with absorption features have gas disks that are seen 
precisely edge-on. 

Keeping this caveat in mind, we assume in the follow- 
ing discussion an outflow origin only for blueshifted HI 
absorption features with Hl-correct ed heliocentric veloc - 
ities (Vq, hJ3) below -50 km s" 1 (|Rupke et al.ll2005aD . 
A number of mechanisms have been proposed to ex- 
plain such HI outflows. They include adiab atically ex- 
pand ing broad emission line clouds (e.g., lElvis et aD 
2002) , starburst- and AGN-driven galacti c winds (e.g., 
Heckman et all UMH I Veilleux et all 120051) . and AGN- 



driven jets interacti ng with the surround ing interstellar 
medium (ISM; e.g.. lOosterloo et aLl l2000). These mech- 
anisms predict different outflow locations (from parsec to 
kiloparsec scales) and geometries (wide angled vs. colli- 
mated). 

6.1. Mass Outflow Rates 
Following iHeckman et al.1 (|2000h and iRupke et al.1 
(2002), the mass outflowing (M) at a velocity v is calcu- 
lated using: 



- 2 ^200 km s" 1 ' 



M yr- 1 1 4tt A 1 kpc A 10 21 cm 

(5) 

assuming spherical geometry, mass conservation, and a 
wind that flows from a distance of to infinity. The out- 
flow is also assumed to subtend a solid angle of Jl within 
which the gas has a covering fraction of 1. To facilitate 
comparison s with the values deriv ed from young radio 
galaxies by IMorganti et al.l (|2005b[ ) , we use the same as- 
sumptions of = 7r and r* = 1 kpc. Applying these as- 
sumptions and Equation 5 to the blueshifted absorption 
components where v = V©, Hi) we find mass outflow rates 

M (wide) up to more than 6000 Mq yr -1 (summarized 
in Table H|). These values are much larger than those 
of <60 M Q yr" 1 measured by IMorganti et~aT1 (|2005bl ) 
and the range me asured in ULIRGs (13-133 M© yr -1 ; 
IRupke et al.ll200l . 

Many of the mass outflow rates we infer seem unreal- 
istically large; such high rates have not been measured 
in other astrophysical objects. These values do, how- 
ever, depend on the assumptions we have made above. 
The mass outflow rate depends linearly on the size scale 
of the outflow region and the derived HI column den- 
sity. Recall from Equation 4 that the column density 

10 These systemic velocities are relative to the HI emission rather 
than the optical emission. See Table [6] for more details. A remain- 
ing source of uncertainty in the interpretation of HI line profiles is 
the question of whether the midpoint of the 20% flux level is always 
a good approximation of the true systemic velocity. Emission from 
molecular lines is a more reliable tracer of systemic velocity than 
HI, so the ideal solution here would be acquisition of single-dish 
CO spectroscopy for the full QUEST sample. Unfortunately, these 
observations are still pending. 



depends on the spin temperature and optical depth pro- 
file, which in turn depends on the brightness of the 
background continuum. There are thus three key fac- 
tors that affect our estimates of a mass outflow rate: 
the spin temperature of the gas, the assumed distance, 
and the background continuum. The assumption of 
1000 K for the HI spin temperature seems to be rea- 
sonably conservative. In the vicinity of powerful AGN, 
the extreme conditions of the gas and the presence of 
shocks can resu lt in spin temperatures >3000 K (e.g., 
Holt et all 120061 ). although we note that absorption-line 
opacity will be dominated by the coldest material in any 
multiphase medium. Spatially resolved interferometric 
maps of several quasars suggest t hat the outflow regions 
have size scales ~0 . 2-1.6 kpc (lOosterloo et al.l 120001: 
Morganti et al.ll2003l: fEmonts et alJl2005UMorganti et al l 



2005aj ). Adopting the smaller distance of 0.2 kpc instead 



of 1 kpc would reduce the derived mass outflow rates 
by a factor of 5, more consistent with those derived for 
ULIRGs. 

Finally, the estimated HI column density also depends 
strongly on the geometry of the source of continuum 
emission. We have assumed that this emission comes 
entirely from the core. However, the absorption troughs 
in several quasars in our sample are deeper than the core 
fluxes. This excess implies that off- nuclear emission also 
contributes to the observed continuum. In this case, t in 
Equation (3) is overestimated and so is Nm, by a factor 
equal to the actual nuclear/total continuum flux ratio, 
although a larger source size may compensate to some 
extent. 

6.2. Jet-induced Outflows 

If the off-nuclear continuum source is indeed physically 
related to the target and not a background quasar that 
happens to lie within the GBT beam, our polarized spec- 
tra give more credence to the idea that radio je ts play an 
impo rtant role in driving HI gas outflows fe.g.. lMorgantil 
l20ll . Several sources (see Appendix) in our sample with 
polarized spectra have spatially resolved radio jets de- 
tected at higher frequencies. The resolutions of these 
VLA data imply that the jets are launched within ~0.4- 
0.9 kpc of the core. Considering a jet-induced outflow 
scenario for the HI absorption features, we again apply 
Equation 5 to the blueshifted outflow components and 
assume r* = 0.5 kpc and a narrow opening angle of 5° 
for the jet. Under these new assumptions, the range 
of mass outflow rates is 0.2-31 Mq yr -1 (see col. (8) 
of Table [5]) . We consider these values lower bounds on 
the actual mass outflow rates, while those derived using 
SI = 7r and ?•» = 1 kpc (col. (7) of Table [6]) are upper 
bounds. 

Under the jet scenario, the kinetic ene rgy release rate 
can be calculated using the equation fromlHcckma n et al.1 
(2000]): 



— ^_ = 2.7xl0 41 ( — )(-£- 
erg s 1 47r 1 kpc 



)(- 



N H i 



200 km s 
(6) 

The kinetic energy in the outflow is on the order of 10 38 
to 10 42 erg s _1 for our sample. 

The association of jets with the detected HI absorption 
features is also seen in the distribution of radio loud- 
ness of our sample. Only one quasar in our sample is 



-)• 



HI in Massive Gas-rich Mergers 



7 



considered radio-loud based on optical-to-radio flux ra- 
tios (see Appendix), but optical data may be subject 
to contamination from the host galaxy, resulting in an 
underestimate of radio loudness. Using the 1.4 GHz to 
2-10 keV ratio as an indicator of radio loudness and 
a radi o-loud cutoff value emp irically derived from 1600 
AGN (|La Franca et al.l 120101) , we find that many ob- 
jects in our sample are radio- loud (Table [T]). These re- 
sul ts are consistent with those found in a recent study 
bv iChandola et ail (|2012l ) that suggest higher detection 
rates of HI absorption toward extended compact steep- 
spectrum and gigahertz-peaked spectrum sources than 
toward the cores of radio galaxies. The environments in 
which the jets interact with the interstellar medium can 
be complex, and we may be seeing a strongly polarized 
continuum being absorbed selectively by multiple neutral 
hydrogen clouds, resulting in polarized HI absorption. 

The fraction of radio-loud sources with associated HI 
absorption (58%) is far higher than the ~10% observed 
in the general population of AGN. Given the high star 
formation rates of the ULIRGs, the circumnuclear star- 
burst can contaminate the 1.4 GHz measuremen ts . Using 
the 1.4 GHz to SFR correlation of Hi3 (|2003t ). we cor- 
rected the 1.4 GHz flux for star formation. Of the three 
ULIRGs with blueshifted absorption and X-ray detec- 
tions, all (100%) are radio-loud. Similarly, of the six 
FIR-strong quasars with detected HI outflows and X-ray 
emission, three (50%) are considered radio-loud; one of 
the two (50%) FIR-weak quasars with HI outflows and 
X-ray flux are radio-loud. This pattern recalls an idea 
put forward by iWilson fc Colbertl JT995) on the forma- 
tion of radio galaxies to explain why radio-quiet AGN 
outnumber radio-loud AGN. The authors suggested that 
major mergers, such as the progenitors of our sample 
of galaxies, result in coalesced supermassive black holes 
with high spin values due to the conservation of angu- 
lar momentum. Substantial energy is then extracted 
from these fast-spinning black holes to drive jets, yielding 
radio-loud AGN. Thus, younger radio galaxies may have 
more powerful jets and therefore stronger jet-induced 
outflow s. This scenario is con sistent with the observa- 
tions of iMorganti et al.l (|2005b| ). whose sample of young 
or recently restarted radio ga l axies s how similar HI out- 
flows. The IWilson fc Colbert (1995) picture is also con- 
sistent with theoretical models of lTchek hovsko v~et al.l 
(2010) that suggest the radio loud/quiet dichotomy in 
AGN is due to the spin of the central black hole. 

6.3. HI Outflows Along the Merger Sequence 

Even though our sample with mass outflow rates is 
relatively small, we see some trends in our outflow mea- 
surements (Table |6]and Figure [12]) . However, given the 
small number of detected HI outflows, particularly with 
only two measurements for the FIR-weak quasars, we 
cannot statically confirm that the distributions we see in 
Figure [T^] is significant. There is also no obvious rela- 
tionships between mass outflow rate and the bolomctric 
luminosity, AGN luminosity, or radio loudness for our 
sample. 

Systems with deep absorption and outflows become 
more prevalent as we progress along the merger sequence 
from the newly coalesced ULIRGs to the FIR-strong 
quasars. The trend between outflow kinematics and star 
formation rates among the ULIRGs (Figures [3] and [TTj) 



suggests that the HI outflows in these objects are driven 
mainly by the starbursts, with a median mass outflow 
rate in the jet scenario (M (jet) ~1 Mq yr" 1 ) similar 
to their SFRs and previous Na I -bas ed outflow measure - 
mcnts (|Rupke et al. l l200l [2005allbl kl: IMartidl200l [2006) . 
ULIRGs also tend to have the narrowest absorption fea- 
tures. In the FIR-strong quasar systems, the median 
mass outflow rate is M (jet) ~5.8 Mq yr -1 , considerably 
larger than their median SFR of ~1.0 Mq yr . While 
the median mass outflow rate of the FIR-weak quasars 
is lower, M (jet) ~0.3 Mq yr -1 , so is their median SFR 
(~0.6 Mq yr -1 ). Thus, the mass-loading factor, i.e., 
the ratio of mass outflow rate to SFR, is much greater 
in the FIR-strong than the value of near unity derived 
in ULIRGs. The same value is much lower in the FIR- 
weak quasars, but this is difficult to compare given that 
this is based on a single measurement. However, the 
distinction between the ULIRGs and FIR-strong quasars 
further supports the idea that the AGN is playing a role 
in driving HI outflows in the quasars before the end of 
large-scale star formation. 

HI in emission is much more prevalent among the 
quasars than the ULIRGs, and the HI emission profiles 
increasingly resemble those of regularly-rotating neu- 
tral gas disks as we move from FIR-strong to FIR- 
weak quasars. The amount of HI seen in emission in- 
creases modestly from ~4xl0 10 Mq among the FIR- 
strong quasars to ^7xl0 10 Mq among the FIR-weak 
quasars. These masses are c onsistent with th ose mea- 
sured in low-redshift quasars (|Konig et al.ll2009fl and the 
high end of the gas m ass range seen in nearby active 
galaxies (|Ho et alil200l . We may be witnessing the cre- 
ation of HI disks from residual neutral gas at the centers 
of merger remnants. 

7. SUMMARY 

We have analyzed of GBT HI spectra of 27 ULIRGs 
and quasars selected from the QUEST sample. We find 
that: 

1. As mergers evolve from fully coalesced ULIRGs 
to FIR-strong quasars to FIR-weak quasars, the 
characteristics of their HI profiles also change from 
mostly absorption in the ULIRGs to mostly emis- 
sion in the FIR-weak quasars, with the FIR-bright 
quasars having intermediate properties. The kine- 
matics of these systems change from narrow ab- 
sorption features likely due to outflows driven by 
starbursts in ULIRGs, to broad absorption features 
from jet-driven outflows driven by AGN in FIR- 
strong quasars, and finally to the emission profiles 
more characteristic of gas in rotation in FIR-weak 
quasars. 

2. A large fraction of the quasars exhibit polarization- 
dependent absorption features, with 78% of FIR- 
strong quasars and 63% of FIR-weak quasars show- 
ing these features compared with 13% in ULIRGs. 
Some of these features vary on short time scales 
(down to ^30 minutes). Neither time variabil- 
ity nor polarization-dependent spectral features are 
seen in contemporaneous observations of M82 (as 
expected), reducing the likelihood that these ef- 
fects are instrumental artifacts. These features are 
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likely due to HI absorption of polarized nuclear and 
off-nuclear continuum emission. The large frac- 
tion (~60%) of radio-loud sources with HI outflow- 
associated absorption suggests a close relationship 
between radio jets and the HI screen. 

3. The polarization-dependent absorption features are 
most likely due to the emergence of radio jets 
in these systems, contributing to the continuum 
emission against which the absorption features are 
detected. The ^60% radio-loud fraction of our 
sources with HI outflows is far higher than the 
~10% seen in the general population of AGN. This 
distinction suggests that radio jets play an impor- 
tant role in shaping the environments in these sys- 
tems, which may tract a transition phase in the 
evolution of gas-rich mergers into "mature" radio 
galaxies. 
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APPENDIX 

NOTES ON INDIVIDUAL GALAXIES 

In this appendix, we describe the HI spectral properties of our GBT sample. When appropriate, we compare the 
present data with previous observations of these galaxies. 

PG 0007+106: The only radio-loud object in our sample based on classic radio-to-optical criteria. This source is 
classified as a triple system in the third Zwicky catalog (III Zw 2), and the observed velocity range is consistent 
with the A (v sys — 26921 km s _1 ) and C (v sys ~ 26981 km s _1 ) components of the system. The GBT detection 
is assumed to be that of the brig htest component, A. Th e HI emission profile of the full intensity spectrum agrees 
with t he old Arecibo spectrum of lHutchings et all (|1987l ). The absorption feature, also noted bv lHutchings et all 
(|1987f ). is more prominent in the XX spectrum of the GBT data. 

PG 0050+124 (I Zw 1): The heliocentric radial velocity listed in Table Q] for this source is from optical spec- 
troscopy bv lHo fc Kim (200jJ- The optically-derived velocity is ^700 km s^ 1 lower than the Hi-derived value 



of 18285 km s" 1 (jBothun et al.|[l~98l . The optical emission lines may arise in part from outflowing gas, so the 
optically-derived value may is underestimate the actual systemic velocity of this obje ct. The shape of th e full 
intensity profile ( Figur e [J) at 18285 km s _1 agrees with those previously published bv lBothun et al. (1984) and 
IHutchings et alj (p.987), and matches expectations for a regularly rotating disk. However, at lower velocities, 
the profile becomes polarized and a broad absorption feature is centered at the optically-derived velocity (see 
Figure©. 

F04103-2838: There is a marginal detection of an emission feature and a weak but broad blueshifted absorption 
feature present. 

F05189-2524: No emission is seen in the data. However, several absorption components are detected. An alternative 
explanation for the spectral shape is the superposition of an emission profile with a broad absorption component. 

PG 0804+761: Strong RFI is present at the expected frequency of this galaxy. It is therefore excluded from our 
analysis. 

PG 0844+249: The broad emission profile is asymmetric. Alth ough the profile s in the two different polarizations 
differ slightly, they have shapes similar to that o bserved bv lHo et al.l (|2008f ). Our integrated GBT flux density 
is larger than that measured by I Ho et all ([2008) . The emission feature redward of the peak may be due to a 
companion at v sys = 19601 km s _1 . 

UGC 05101: Multiple absorption components are seen in the GBT data. We interpret the HI absorption profile as 
a highly absorbed core coincident with an inclined regularly rotating disk. 

PG 1119+120: This object was previously observed by iBothun et ail (|1984D and IHutchings et all (|1987l ) with 
Arecibo. The profiles of the older spectra resemble our GBT observation, although the narrow absorption 
feature is not apparent in the lower-resolution Arecibo data. There is n o obvious detection of a companion in the 
GBT data near 16000 km s _1 as suggested by IHutchings et al.l (|1987l ). High resolution VIA maps at 4.8 GHz 
show a bent jet extending to the northeast (|Leipski et al.ll2006l ). At the resolution of these VLA maps (^0.75"), 
the jet is resolved at a distance of ~0.72 kpc. 

PG 1126-041: This target was statistically undetected after three hours on-source, although there is a very weak 
broad absorption feature centered near its systemic velocity. 

PG 1211+143: Thi s broad absorption line system was undetected bv lCondon et alj (|1985f ); when it was observed by 
IHo et al.1 (|2008l ) with Arecibo, standing waves were present in the data due to the strong continuum. 

PG 1229+204: This object was observed bv lHo et all (|2008D . The profile of the GBT spectrum roughly agrees with 
the earlier data. 

PG 1244+026: A broad, blueshifted absorption feature is present. As for F05189-2524, an alternative explanation 
for the spectral shape may be emission absorbed by a single broad absorption feature. 
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Mrk 231: The deep abs orption trough is associated with the face-on disk of the galaxy (|Carilli et al.l 119981 : 
lUlvestad et "all 1 19991 ) . Our GBT full intensity HI profile of this sources mostly agrees with these previous 
observations, but we detect an a dditio nal absorption feature centered at relative velocity ~-300 km s _1 that 
was first suspected by iMorgantil (| 2 1 If ) . This feature is slightly polarized, with the XX spectrum having lower 
intensity. Although the polarization is not as high as those seen in Figure [4] if we give the spectra of Mrk 231 
the same treatment detailed in Section [5721 we find that this feature has FWHM ^1000 km s _1 corresponding 
to a mass out flow rate of ~0-3 yr" 1 in the jet sc enario. The polarized feature is thus likely associated with 
the radio jet (jCarilli et al.lll998l : lUlvestad et al.lfl999f ) . The velocity of this jet-enhanced outflow is similar to the 
value measured using Na I in the jet-disturbed region (-1400 km s-^ IRupke fc VeilleuxllMll) . 

Mrk 273: This sour ce is the only binary in our sample; its nuclear separation is ~750 kpc (jScoville et al.|[2000l) . X-ray 
observations by|Twasawa ct al. (201l|) identified its southwestern nucleus as the Seyfert 2 nucleus and suggested 
that the northern nucleus may host a heavily obscured AGN. The double-peaked HI absorption profile is very 
symmetric and is likely due to the two obscured nuclei. There appears to be a weak absorption feature in the 
YY spectrum at -500 km s _1 . 

PG 1351+640: The XX spectr um appears to show emission, while YY shows broad absorption. The HI absorption 
data may also be variable. iLeipski et al.l (|2006l ) note a jet -like structure in th eir 4.8 GHz VLA maps and time 
variability in the core flux density. The resolution of the ILeipski et al.l ()2006l ) data (~0.56") corresponds to a 
distance of ^0.91 kpc. Although on a different spatial scale, this structure is consistent (si milar position angle) 
with a single-sided extension hinted at by the 8.4 GHz VLB A snapshot data presented bv iBlundell fc Beaslevl 
(1998). This feature is not seen in deeper VLBA imaging at multiple frequencies, but it should be noted that 
the data from lUlvestad et al.l (2005) are only in the left-circular polarization. 

PG 1411+442: The spectrum of this object is very similar to that of PG 1351+640. Emission is seen in the XX 
spectrum and a two-component absorption feature in the YY spectrum. 

PG 1426+015: iHo et al.l (|2008D detect a fairly narrow (W 20 = 357.4 km s" 1 ) and well-defined HI line centered at 
v sys = 25975.2 kms -1 . We have examined our spectrum at the corresponding frequency and see hints of emission 
over a wider range of velocities; however, this feature has such low significance and is so sensitive to the details 
of baseline fitting we classify its detection status as "uncertain." 

PG 1440+356: The profile of the HI absorption feature depends slightly on the polarization. Exceptionally, the 
feature is redshifted by ^300 km s^ 1 with respect to systemic velocity rather than being blucshifted. A possible 
explanation is that the systemic velocity taken from NED and based on optical emission lines underestimates the 
actual systemic velocity because the line-emitting gas arises in part from outflowing material. Another possibility 
is that we may be detecting infalling HI. Emission is also present in the blue half of the spectrum, although this 
is only seen at one of several epochs. 

PG 1448+273: The weakly polarized spectra show blueshifted absorption. 

PG 1501+106: A weak emission feature is present in the data. A blueshifted absorption feature is also visible, 
particularly in the XX spectrum. A lot of the data for this source were affected by small ripples in the baseline. 
As a reality check, we compared our data with similar GBT data from program 10A-070 (PI: Kannappan). The 
ripples also appear to be present in a large fraction of their data. The persistent ripples could be due to harmonic 
waves from a strong RFI signal at a frequency that is slightly outside the detector frequency range. Due to the 
expected weakness of the features, we conservatively flagged all of the data that showed these ripples. 

F15130-1958: Multiple weak absorption components are present in this object. 

F15250+3608: This is the only ULIRG in our sample that shows distinct polarization-dependent HI absorption. The 
spectra are very similar to those of PG 1351+640 and PG 1411+442, where weak emission is seen in the XX 
spectrum and absorption in YY. 

F15462-0450: Strong RFI is present at the expected observed frequency of this galaxy. It is therefore excluded from 
our analysis. 

PG 1617+175: Broad emission is detected in this object. The profiles of the emission differ slightly in the two 
polarizations. Notably, the emission feature in the YY spectrum seems to be affected by an additional absorption 
feature (Figure [5]). 

F21219-1757: This object was not detected. 

PG 2130+099: This target was observed bv lHo et al.l (|2008l ). Our full intensity profile of this s ource agrees w i th the 
older data from Arecibo. A blueshifted absorption feature is also seen in the XX spectrum. ILeipski eF al. (2006) 
detected what looks like jet-associated hot spots emanating from the core in 4.8 GHz A-array VLA data. The 
resolution of these VLA data (^0.35") implies that the jet can be resolved at a distance of ^0.42 kpc. 

PG 2214+139: The HI spectrum of this object is dominated by a broad emission feature centered at systemic velocity. 



TABLE 1 

The Sample and Journal of Observations 



Source 


1 , A 

HA 


Dec 


f-'sys 


Type 


FIR 


l og (_tei) 


SFR 


Si. 4 GHz 


log Rx 


Dates of 


tint 


AT/ 




Fit 




J2000 


J2000 


(km s — 1 ) 




Strength 


(M Q yr" 1 ) 


(Jy) 




Observations 


00 


(km s" 1 ) 


(mjy) 


Order 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


PG 0007+106 


00:10:31.0 


+10:58:30 


26783 


QSO 


Weak 


12.23 


~1.2 


0.099 


-2.4 


2011 Sept - 2012 Jan 


7894 


6.11 


0.77 


3 


FG 0050+124 


00:53:34.9 


i 1 o At o r 1 

+ 12:41:36 


17658 


QSO 


Strong 


12.07 


10.4 


0.009 1 


-5.0 


OOI 1 CI j_ 001 O T 

2011 Sept - 2012 Jan 


4815 


5.78 


0.94 


3 


F04103-2838 


04:12:19.4 


-28:30:25 


35215 


ULIRG 


Strong 


12.30 


154.2 


0.013 1 


-3.4 


2011 Sept - 2012 Mar 


7452 


6.43 


0.88 


3 


F05189-2524 


05:21:01.4 


-25:21:45 


12760 


ULIRG 


Strong 


12.22 


71.6 


0.029 1 




2011 Sept - 2012 Jan 


12330 


5.60 


0.58 




PG 0804+761 


08:10:58.6 


+76:02:43 


29979 


QSO 


Strong 


12.08 


0.3 


0.004 1 


-5.3 


2011 Aug 


860 


6.23 


2.97 




PG 0844+ 349 


o o A *~7 /i O A 

08:47:42.4 


1 o A A P" O A 

+34:45:04 


19187 


QSO 


Weak 


11.44 


0.8 


^ r\ 000 r" 1 

< 0.0025 1 


<-3.3 


2011 Nov 


3153 


5.83 


0.96 


3 




Uy.oo.ol.u 


i ci .01 .1 1 
+ol.zi.ll 


1 1 Qno 
lloUz 


TTT TT3P 


Strong 


1 nc; 
IZ.Uo 


linn 
11U.2 


n 1 7n 
U. 1 / U 


q n 
— o.U 


zUll Aug — zUll ocpt 


Ol c: -| 
Zlol 


o.o I 


1 AP. 

1.40 


Q 
O 


r\j 1119+12U 


11 1 /I7 1 

11:21:47.1 


I 1 1 A A . 1 O 

+ 11:44:18 


15050 


QSO 


Strong 


11.33 


0.4 


0.004 1 




2011 Oct — 2012 Mar 


3837 


5.68 


0.83 


3 


PG 1126-041 


11:29:16.6 


-04:24:08 


18575 


QSO 


Strong 


11.52 


2.5 


< 0.0025 1 


<-4.6 


2011 Aug - 2012 Jan 


8871 


5.81 


0.54 


3 


Y\j lzll+14o 


1(1.1 ,1.17 i-7 

lz:14:l { . / 


1 1H.HQ.1Q 

+ 14:Uo:lo 


24253 


QoO 


Strong 


11.96 


0.0 


0.600 


—2.8 


zOll Aug — zUlz Mar 


5986 


6.02 


0.92 


3 


dp 1 oon i on/i 
Jt^Lj lZZy+ZU4 


1 O- QO.DQ A 

IZ.oZ.Uo.u 


+zu.uy.zy 


looyU 




\\ oak 


11.00 


n A 
{). 1 


n nn qI 
U.UUo 


1 
— 4.y 


ZUlz ivlar 


o411 


qo 
O.oz 


n no 


Q 
O 


r^Lr 1Z44+UZO 


1Z.40.O0.Z 


+uz. zz. uy 


1 A A At 
1444o 




Strong 


1 1 no 
11. (JZ 


1.0 


n nno2 
U.UUz 


—0.1 


zuii Aug — zuiz ivlar 


Qt; 

OOOO 


O.OD 


l.UU 


Q 

o 




1 9-^FM A 9 
1Z . OU. -L^i. Z 


-|-ou. oz.zo 


12642 


TTT TRP 


Strong 


1 9 fin 


174 2 


u.ouy 


—9 7 

Z. 1 


901 1 Alio- 901 1 ^ont 


1 9^7 
1ZO / 


^ RO 
o.uu 


1 52 


Q 
O 


Mrk 273 


13:44:42.1 


+55:53:13 


11326 


ULIRG 


Strong 


12.24 


141.6 


0.136 


-3.5 


2011 Aug - 2011 Sept 


1290 


5.50 


1.31 


3 


PG 1351+640 


13:53:15.8 


+63:45:46 


26442 


QSO 


Strong 


12.04 


26.1 


0.034 1 


-3.3 


2011 Sept - 2012 Feb 


8153 


6.10 


0.79 


3 


PG 1411+442 


14:13:48.3 


+44:00:14 


26861 


QSO 


Weak 


11.78 


0.0 


< 0.0025 1 


<-4.6 


2011 Aug - 2012 Feb 


6346 


6.12 


0.92 


3 


PG 1426+015 


14:29:06.6 


+01:17:06 


25953 


QSO 


Weak 


11.92 


1.9 


0.004 1 


-5.0 


2012 Mar 


3153 


6.08 


1.09 


5 


PG 1440+356 


14:42:07.4 


+35:26:23 


23700 


QSO 


Strong 


11.80 


10.3 


0.005 1 


-4.8 


2011 Aug - 2011 Dec 


6135 


6.00 


0.63 


3 


PG 1448+273 


14:51:08.7 


+27:09:27 


19487 


QSO 


Weak 


11.43 


0.1 


0.003 3 


-4.8 


2011 Aug - 2011 Feb 


8917 


5.85 


0.63 


3 


PG 1501+106 


15:04:01.2 


+10:26:16 


10919 


QSO 


Strong 


11.33 


0.0 


< 0.0025 1 


<-5.7 


2012 Feb - 2012 Mar 


6540 


5.53 


0.77 


3 


F15130-1958 


15:15:55.2 


-20:09:17 


32585 


ULIRG 


Strong 


12.23 


108.8 


0.010 1 


-2.4 


2011 Dec - 2012 Mar 


2666 


6.33 


1.25 


3 


F15250+3608 


15:26:59.4 


+35:58:38 


16535 


ULIRG 


Strong 


12.12 


109.6 


0.018 




2011 Aug - 2011 Nov 


2007 


5.74 


1.16 


. 3 . ( 


F15462-0450 


15:28:56.8 


-04:59:34 


29917 


ULIRG 


Strong 


12.28 


168.7 


0.013 1 


-3.4 


2011 Dec 


860 


6.23 


4.14 




PG 1617+175 


16:20:11.3 


+17:24:28 


33708 


QSO 


Weak 


11.74 


4.0 


< 0.0025 1 




2011 Aug - 2012 Mar 


2547 


6.38 


1.59 


3 


F21219-1757 


21:24:41.6 


-17:44:46 


33526 


ULIRG 


Strong 


12.17 


48.3 


0.015 1 




2011 Oct - 2012 Mar 


8529 


6.37 


0.70 


3 


PG 2130+099 


21:32:27.8 


+10:08:19 


18880 


QSO 


Strong 


11.77 


0.2 


0.007 1 


-2.7 


2012 Mar 


2001 


5.82 


1.13 


5 


PG 2214+139 


22:17:12.2 


+14:14:21 


19715 


QSO 


Weak 


11.77 


0.0 


< 0.0025 1 


<-5.0 


2011 Sept - 2011 Oct 


2121 


5.85 


1.32 


3 



References. — ^Condon etU] fNVSS; fT99l ); HRafter et "all 120091) ; IBarvainis et al.| fl99l) . 

Note. — Col.(l): Galaxy name. Coordinate-based names beginning with "F" are sources in the IRAS Faint Source Catalog; the "PG" sources are Palomar Green quasars. Col. (2): Right 
ascension in J2000 coordinates. Col. (3): Declination in J2000 coordinates. Col. (4): Heliocentric radial velocity given by NED. These velocities arc measured via optical emission lines and may 
underestimate the true recessional velocities since some fraction of t he line emission ma y arise from the outflowing material itself. All velocities are quoted in the optical convention where v — cz. 
Col. (5): ULIRG or QSO source. Col. (6): FIR strength as defined by Nctzcr et al. (2007) where FIR-strong sources have /eo^m/Zis^m > 1- Col. (7): Bolometric luminosity. For ULIRGs, we assume 
L bol — 1.15-Lir. For PG QSOs, we assume L^ol — 7L(5100 A) (Nctzcr et al. 2007). Col. (8): Star formation rates in solar masses per year derived using the average AGN contribution fraction to 
the total bolometric luminosity measured by IVeilleux et al.l d2009al) for each galaxy. For PG 0007+106, the AGN fraction is assumed to be the me an for all o t her FIR-wcak quasars. It is assumed 
that all of the emission from the starburst is reprocessed by dust and re-emitted into the IR band. The SFR is derived using the relation given in IKennicutt] U998I) for a 0.1-100 M initial in ass 
function. Col. (9): 1.4 GHz continuum flux densities. Values taken with the GBT unl ess noted. The upper limits are given as the detection limit of the NRAO VLA Sky Survey (NVSS). Col. (10): 
1.4 GHz to 2—10 k eV flux ratio, where lo g Rx > —4.3, implies a radio-loud source I ILa Franca et al.| [2010). We used the (avera ge, if multiple observa tions were performed) absorption-corrected 
2-10 keV flux from lTeng & Veilleuxl d2010fl except for PG 0007+106 and PG 0804+761 whose X-ray meas urements w ere taken from lPiconcelli et all d2005li . These ratios are correct ed for the 14 GH z 
contribution from star formation based on the SFRs listed in Tableland the SFR-1.4 GHz relation from Bell (2003), which assumes the same 0.1—100 Mq initial mass function as Kcnnicutt ( 1998). 
For ULIRGs, due to their high line-of-sight column densities, their measured 2-10 keV fluxes may be underestimated and thus the listed values may represent upper limits. Col. (11): Period over 
which the target was observed. Col. (12): Total on source time in seconds after flagging. Col. (13): Velocity resolution of the full intensity spectrum. Col. (14): The channcl-to-channel rms of the full 
intensity spectrum. Col. (15): The order of polynomial over which the baseline was fitted and then subtracted. 
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TABLE 2 

Measurements of HI in Emission 



Source 


y@ 


W20 


w 50 


Peak 


S/N 


/hi 




N m 


M m 




(km s _1 ) 


(km s — 1 ) 


(km s" 1 ) 


(mjy) 




(Jy km s 


_1 ) 


(10 18 cm" 2 ) 


(10 10 Mq) 


(1) 


(2) 


(3) 


(4) 


(5) 


((>) 


(7) 




(8) 


(9) 


PG 0007+106 


-216 + 16 


1063 ± 49 


894 ± 33 


4.30 


5.58 


2.57 + 0. 


.14 


9.35 ± 0.49 


8.46 ± 0.45 


PG 0050+123 


473 ± 17 


789 ± 51 


435 ± 34 


7.00 


7.45 


2.43 + 0. 


14 


8.85 + 0.51 


3.44 + 0.20 


F04103-2838 (M) 


129 + 21 


151 ± 63 


21 + 42 


3.44 


3.91 


0.17 + 0. 


06 


0.62 ± 0.22 


0.97 + 0.34 


PG 0844+249 


243 + 8 


541 ± 23 


436 ± 15 


8.83 


9.20 


2.32 + 0. 


12 


8.44 ± 0.43 


4.05 + 0.21 


PG 1119+120 


-18 + 2 


226 ± 7 


221 + 5 


5.35 


6.45 


0.58 + 0. 


07 


2.11 ± 0.24 


0.65 ± 0.07 


PG 1211+143 


522 + 8 


200 ± 25 


174 ± 17 


3.85 


4.18 


0.04 + 0. 


07 


0.14 + 0.25 


0.11 + 0.20 


PG 1229+204 


182 ± 11 


600 ± 34 


504 ± 23 


5.37 


5.84 


1.84 + 0. 


12 


6.70 ± 0.43 


3.24 + 0.21 


Mrk 273 


407 ± 11 


116 ± 32 


62 + 21 


5.94 


4.53 


0.37 + 0. 


07 


1.35 + 0.26 


0.24 + 0.05 


PG 1351+640 


-205 ± 21 


585 ± 64 


480 ± 43 


2.64 


3.34 


0.82 + 0. 


10 


2.98 ± 0.38 


2.73 ± 0.34 


PG 1426+015 (U) 


107 ± 35 


1131 ± 104 


591 ± 69 


5.12 


4.69 


2.08 + 0. 


20 


7.57 + 0.72 


6.86 + 0.65 


PG 1440+356 


-294 ± 13 


477 ± 39 


395 ± 26 


3.01 


4.81 


0.48 + 0. 


07 


1.73 ± 0.27 


1.30 + 0.20 


PG 1448+273 


-50 ± 18 


97 + 54 


23 + 36 


2.04 


3.26 


0.11 + 0. 


03 


0.41 + 0.12 


0.21 + 0.06 


PG 1501+106 


61 ± 21 


137 ± 63 


22 + 42 


2.58 


3.37 


0.32 + 0. 


05 


1.16 + 0.17 


0.20 + 0.03 


PG 1617+175 (U) 


-118 + 15 


1139 + 44 


1057 ± 30 


6.94 


4.36 


3.88 + 0. 


30 


14.12 ± 1.08 


21.21 ± 1.62 


PG 2130+099 


224 + 9 


599 ± 26 


518 ± 18 


7.86 


6.95 


2.30 + 0. 


15 


8.37 + 0.53 


3.84 + 0.24 


PG 2214+139 


104 ± 14 


876 ± 42 


567 ± 28 


11.22 


8.53 


4.43 + 0. 


21 


16.13 + 0.75 


7.99 + 0.37 



Note. — Col.(l): Galaxy name where (M) and (U) indicate marginal and uncertain detections, respectively, as noted in 
Section |3.2| Col. (2): Heliocentric velocity measured from HI relative to the systematic velocity presented in Table HI and 
its associated measurement error. Col. (3): Velocity width of the emission measured at the 20% level and its error. Col. (4): 
Velocity width of the emission measured at the 50% level and its error. Col. (5): Maximum intensity of the spectrum. Col. (6): 
Signal-to-noise ratio given by the peak intensity relative to the channel rms listed in Table [T] Col. (7): Velocity-integrated HI 
flux in emission and its error. Col. (8): Derived HI column density calculated using Equation 1. Col. (9): Derived HI gas mass 
from emission and its error calculated using Equation 2. 
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TABLE 3 

Measurements of HI in Absorption 



Source 


^3 


FWHM 


"Peak" 


S/N 




Tin ax 


Jabs 




(km s _1 ) 


(km s _1 ) 


(mjy) 


(10 : 


22 cm ~ 2 ) 




(Jy km s _1 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 






Strong Detections: peak-to 


-rms > 3er 








F05189-2524 


160 + 3 


79 + 8 


-2.74 + 0.23 


4.73 


1.50 


0.099 


-0.04 


UGC 05101 


-251 ± 11 


419 ± 26 


-6.69 + 0.22 


4.59 


3.24 


0.040 


-0.54 




8 + 3 


96 + 7 


-9.76 + 0.59 


6.70 


1.09 


0.059 


-0.18 




203 + 9 


219 ± 22 


-5.23 + 0.34 


3.59 


1.32 


0.031 


-0.22 


PG 1119+120 


-88 + 3 


22 + 6 


-4.20 + 0.87 


5.06 






-0.02 


PG 1211+143 


-227 + 1 


460 + 3 


-37.67 + 0.16 


40.95 


5.73 


0.065 


-3.07 


PG 1244+026 


-381 ± 8 


226 ± 20 


-3.07 + 0.23 


3.05 






-0.13 


Mrk 231 


-6 + 2 


199 + 6 


-16.87 + 0.41 


11.07 


2.15 


0.056 


-0.64 


Mrk 273 


-122 + 8 


334 ± 21 


-11.41 + 0.41 


8.71 


5.60 


0.088 


-0.74 




187 + 6 


171 ± 13 


-10.62 + 0.64 


8.11 


2.66 


0.081 


-0.35 


PG 1440+356 


325 + 4 


333 + 8 


-9.57 + 0.21 


15.27 






-0.57 






Weak Detections: peak-to- 


rms < 3cr 








nr 1 nnc;n i io/l 
r\j UUdU+124 


15 ± 10 


53 + 24 


-1.34 + 0.51 


1.42 


1.62 


0.161 


—0.01 


T?C\A 1 HQ /T\ /T \ 

r U41Uo— zooo (^JVl J 


i no _L on 
—409 ± 29 


449 ± 70 


-1.17 + 0.16 


1.33 


8.10 


0.094 


—0.09 


T?nc;i on oko/I 
r Uoloy— zbz4 


—374 ± 9 


118 + 22 


-1.29 + 0.20 


2.22 


1.03 


0.046 


—0.03 




—79 ± 13 


163 ± 30 


-1.09 + 0.17 


1.88 


1.20 


0.038 


-0.03 




347 ± 11 


127 ± 27 


-1.05 + 0.19 


1.81 


0.90 


0.037 


-0.03 


"DO Inn i 1 /I o 


a a a _l_ o e; 
— 644 ± 20 


214 ± 60 


-1.14 + 0.25 


1.24 


0.08 


0.002 


—0.04 




766 + 11 


116 + 26 


-2.47 + 0.48 


2.68 


0.09 


0.004 


—0.05 


PG 1244+026 


183 + 18 


268 ± 43 


-1.52 + 0.21 


1.52 


59.07 


1.407 


—0.08 


Mrk zol 


—280 + 9 


96 + 21 


-3.14 + 0.58 


2.06 


0.19 


0.010 


—0.06 


rij 1351+640 


nil a _1_ 1 n 

— ZZ4 ± 10 


115 ± 25 


-1.53 + 0.28 


1.94 


1.02 


0.046 


—0.03 




336 ±61 


516 ± 98 


-0.97 + 0.14 


1.23 


2.88 


0.029 


—0.09 


PG 1411+442 


-446 ± 10 


319 ± 23 


-2.60 + 0.16 


2.81 






-0.14 




-41 ± 10 


214 + 25 


-1.70 + 0.17 


1.84 


>39.78 


1.139 


-0.06 




319 + 10 


143 ± 24 


-1.76 + 0.26 


1.90 


>28.08 


1.217 


-0.04 


PG 1448+273 


-441 ± 14 


392 ± 34 


-1.69 + 0.13 


2.70 


55.58 


0.829 


-0.12 




144 ± 19 


195 ± 45 


-0.74 + 0.15 


1.18 


10.27 


0.283 


-0.03 


PG 1501+106 


-229 + 12 


304 ± 28 


-2.12 + 0.17 


2.77 


>83.12 


1.884 


-0.12 




340 ± 20 


195 ± 48 


-0.98 + 0.21 


1.28 


>17.46 


0.498 


-0.04 


F15130-1958 (M) 


-847 ± 10 


89 + 23 


-2.07 + 0.47 


1.66 


3.86 


0.232 


-0.03 


-451 ± 15 


168 ± 36 


-1.75 + 0.32 


1.40 


6.10 


0.192 


-0.05 




199 ± 16 


182 ± 38 


-1.74 + 0.32 


1.40 


6.56 


0.191 


-0.05 


F15250+3608 


-42 ± 8 


83 ± 18 


-2.57 + 0.49 


2.21 


2.42 


0.154 


-0.04 




229 + 8 


107+19 


-2.79 + 0.43 


2.40 


3.41 


0.168 


-0.06 



Note. — Col.(l): Galaxy name. Col. (2): Heliocentric velocity measured from HI relative to the systematic 
velocity presented in Table [T1 as modeled using a Gaussian model, and its associated measurement error. Col. (3): 
Full width of the absorption feature at the 50% flux level as measured using a Gaussian model and its error. 
Col. (4): Peak absorption flux and its error as measured using a Gaussian model. Col. (5): Signal-to-noise ratio 
measured as a ratio between the peak flux and channel rms listed in Tabled Col. (6): Integrated HI column density 
based on the absorption model, assuming a spin temperature of 1000 K and using Equation 4 based on the optical 
depth profile derived using Equation 3. For PG 1411+442 and PG 1501 + 106, the 1.4 GHz core flux densities are 
only upper limits, so their derived -/Vhi values are lower limits. Col. (7): Peak optical depth based on Equation 3. 
Col. (8): Velocity-integrated HI flux in absorption measured from the best-fit Gaussian component models. 
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TABLE 4 

Measurements of HI in Emission from Complex Spectral Modeling 



Source 


y@ 


W 20 


w 50 


"Peak" 


S/N 




fin 




N m 


M m 




(km s — 1 ) 


(km s — 


(km s — ) 


(mJy) 




(Jy km s 


_1 ) 


(10 18 cm" 2 ) 


(10 10 M G} ) 


(!) 


(2) 


(3) 


(4) 


(5) 


(6) 




(7) 




(8) 


(9) 


PG 0007+106 


-211 ± 11 


1054 ± 33 


917 + 22 


5.74 


7.41 


3. 


81 ± 0. 


.14 


13.86 + 0.50 


12.53 ± 0.45 


PG 0050+123 


184 ± 21 


1367 ± 64 


916 ± 42 


8.32 


6.80 


i. 


53 ± 0. 


.24 


16.48 + 0.87 


6.40 ± 0.34 


PG 0844+249 


120 ± 14 


921 ± 43 


474 ± 29 


10.65 


10.06 


1. 


.36 ± 0. 


17 


15.86 + 0.62 


7.61 ± 0.30 


PG 1119+120 


-15 + 3 


232 + 8 


222 + 6 


6.23 


7.43 


0. 


.90 ± 0. 


07 


3.28 + 0.24 


1.02 ± 0.08 


PG 1211+143 


-235 ± 17 


1802 ± 50 


655 ± 34 


13.10 


13.97 


7. 


.89 ± 0. 


.21 


28.73 ± 0.78 


22.61 + 0.67 


PG 1351+640 


-332 ± 29 


1202 ± 88 


807 ± 59 


3.56 


4.72 


1. 


.34 ± 0. 


.14 


4.89 + 0.52 


4.48 ± 0.47 


PG 1411+442 


-205 ± 14 


208 ± 41 


169 ± 27 


3.22 


3.19 


0. 


.42 ± 0. 


OS 


1.52 + 0.29 


1.45 + 0.28 


PG 1448+273 


-68 ± 12 


168 ± 37 


91 ± 25 


3.07 


4.83 


0. 


.26 + 0. 


04 


0.93 + 0.16 


0.48 ± 0.08 


PG 1501+106 


61 + 5 


135 ± 16 


122 ± 11 


3.22 


4.53 


0. 


.22 + 0. 


04 


0.81 + 0.16 


0.14 + 0.03 


F15250+3608 


-89 ± 10 


1151 + 30 


1122 ± 20 


5.44 


3.67 


2. 


.38 ± 0. 


.26 


8.66 + 0.96 


3.23 + 0.36 


PG 1617+175 (U) 


-88 ± 13 


1315 ± 38 


1208 ± 25 


9.23 


5.81 


6. 


.10 + 0. 


.32 


22.19 + 1.16 


33.32 ± 1.75 


PG 2130+099 


67 + 20 


909 ± 61 


366 ± 40 


8.56 


7.87 


2. 


.62 ± 0. 


17 


9.54 + 0.63 


4.37 + 0.29 



Note. — Col.(l); Galaxy name. Col. (2): Heliocentric velocity measured from HI relative to the systematic velocity presented 
in Tabled and its associated measurement error. Col. (3): Velocity width of the emission measured at the 20% level and its 
error. Col. (4): Velocity width of the emission measured at the 50% level and its error. Col. (5): The maximum intensity of 
the spectrum. Col. (6): Signal-to-noisc ratio given by the peak intensity relative to the channel rms listed in Table [T1 Col. (7): 
Velocity-integrated HI flux in emission and its error. Col. (8): Derived HI column density calculated using Equation 1. Col. (9): 
Derived HI gas mass from emission and its error, calculated using Equation 2. 



TABLE 5 

Measurements of HI in Absorption from Complex Spectral Modeling 



Source 


v© 


FWHM 


"Peak" 


S/N 


Nm 


7"max 


/abs 




(km s 1 ) 


(km s" 1 ) 


(mJy) 




(10 22 cm ~ 2 ) 




(Jy km s _1 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


PG 0007+106 


-270 ± 18 


668 ± 43 


-3.26 + 0.18 


3.28 


4.32 


0.033 


-0.38 


PG 0050+124 


-32 ± 22 


894 ± 58 


-4.63 + 0.23 


3.41 


112.44 


0.722 


-0.76 


PG 0844+249 


158 ± 22 


700 ± 51 


-4.52 + 0.29 


2.42 






-0.58 


PG 1119+120 


-82 + 2 


25 + 4 


-8.49 + 1.25 


4.73 






-0.04 




17 + 27 


244 ± 56 


-2.00 + 0.40 


1.11 


29.66 


0.693 


-0.09 


PG 1211+143 


-270 ± 4 


466 ± 14 


-37.80 ± 1.88 


20.43 


5.82 


0.065 


-3.11 




-61 ± 12 


870 ± 12 


-36.10 + 1.77 


19.52 


10.38 


0.062 


-5.56 


PG 1351+640 


-68 ± 18 


952 ± 45 


-3.77 + 0.15 


3.13 


21.32 


0.118 


-0.63 


PG 1411+442 


-288 + 11 


735 ± 27 


-6.19 + 0.20 


5.21 






-0.79 


PG 1448+273 


-106 + 29 


723 ± 69 


-1.98 + 0.16 


1.54 


128.80 


1.079 


-0.26 


PG 1501+106 


-140 + 27 


550 ± 64 


-2.15 + 0.22 


1.47 


>155.41 


1.966 


-0.23 


F15250+3608 


-71 ± 21 


858 ± 52 


-5.76 + 0.29 


2.73 


60.56 


0.386 


-0.92 


PG 1617+175 (U) 


-95 ± 30 


870 ± 71 


-4.60 + 0.32 


2.08 






-0.67 


PG 2130+099 


187 ± 35 


325 ± 87 


-2.52 + 0.55 


0.83 


26.35 


0.446 


-0.15 




-247 ± 14 


172 ± 33 


-4.68 + 0.74 


1.55 


31.27 


1.104 


-0.15 



NOTE. — Col.(l): Galaxy name. Col. (2): Heliocentric velocity measured from HI relative to the systematic 
velocity presented in Tabled as modeled using a Gaussian fit, and its associated measurement error. Col. (3): Full 
width of the absorption feature at the 50% flux level as measured using a Gaussian model and its error. Col. (4): 
Peak absorption fiux and its error as measured using a Gaussian model. Col. (5): Signal-to-noisc ratio measured as 
a ratio between the peak flux and channel rms listed in Table Col. (6): Integrated HI column density based on 
the absorption model, assuming a spin temperature of 1000 K, and using Equation 4 based on the optical depth 
profile derived using Equation 3. For PG 1501 + 106, the 1.4 GHz core flux density is only an upper limit, so its 
derived -/Vhi value is a lower limit. Col. (7): Peak optical depth based on Equation 3. Col. (8): Velocity-integrated 
HI flux in absorption measured from the best-fit Gaussian component models. 
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TABLE 6 

Measurements of HI Outflow 



Source 


V & 


V©, HI 


Nm 


FWZI 


Vmax 


^98% 


M (wide) 


M (jet) 


, / B(jct) n 




(km s _1 ) 


(km s _1 ) 


(10 21 cm " 2 ) 


(km s" 1 ) 


(km s" 1 ) 


(km s" 1 ) 


(Af Q yr" 1 ) 


(M yr" 1 ) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


ULIRGs 


F04103-2838 (M) 


-409 


-538 ± 36 


81 


1267 + 212 


-1171 ± 214 


-920 ± 70 


1089 


4.4 


41.6 


F05189-2524 


-374 


-374 ± 9 


10 


336 ± 67 


-542 ± 68 


-475 ± 21 


97 


0.4 


40.3 




-79 


—79 ± 13 


12 


459 ± 90 


—309 ± 91 


—218 ± 29 


24 


0.1 


38.3 
















121+ 


0.5t 




Mrk 231 


-280 


-280 ± 9 


2 


297 ± 73 


-428 ± 74 


-361 ± 20 


13 


0.1 


39.1 


F15130-1958 (M) 


-847 


-847 ± 10 


39 


266 ± 76 


-980 ± 77 


-923 ± 22 


818 


3.3 


41.9 


-451 


-451 ± 15 


61 


494 + 114 


-697+ 115 


-594 ± 34 


687 


2.7 


41.3 
















1505+ 


6.0t 




F15250+3608 


-71 


18 ± 23 


606 


2744 ± 172 


-1354 ± 173 


-53 ± 28 










Median 




39 


459 ± 90 


-697 ± 91 


-475 ± 28 


605 


1.6 


40.8 










FIR-strong 


Quasars 










PG 0050+124 


-32 


-216 + 30 


1124 


2815 ± 191 


-1623 + 192 


-977 ± 58 


6072 


24.3 


41.6 


PG 1119+120 


-82 


-67 + 4 




85 ± 17 


-110 + 17 


-89 + 5 










17 


32 + 27 


297 


721 ± 182 


-329 ± 184 


-176 + 55 








PG 1211+143 


-270 


-270 ± 17 


58 


1680 ± 54 


-1109 + 54 


-666 ± 21 


392 


1.6 


40.6 




-61 


-61 ± 21 


101 


3118 ± 55 


-1620 + 56 


-801 ± 23 


158 


0.6 


38.9 
















550 1 " 


2.2 1 " 




PG 1244+026 


-381 


-381 ± 18 


152 


696 ± 62 


-729 ± 63 


-573 ± 25 


1450 


5.8 


41.4 


PG 1351+640 


-68 


264 ± 34 


213 


2946 ± 147 


-1210 + 148 


-547 ± 51 








PG 1440+356 


325 


325 + 4 




1104 + 45 


-227 ± 45 


41 + 8 








PG 1501+106 


-140 


-201 ± 27 


> 1554 


1637 ± 199 


-1020 + 201 


-669 ± 61 


> 7813 


> 31.3 


> 41.6 


PG 2130+099 


-247 


-314 ± 24 


313 


547 + 111 


-587+ 112 


-460 ± 37 


2453 


9.8 


41.5 




Median 




213 


1370 ± 87 


-1020 + 88 


-560 ± 31 


1450 


5.8 


41.5 










FIR-weak 


Quasars 










PG 0007+106 


-270 


-59 + 21 


43 


2053 ± 133 


-1085 ± 134 


-627 + 42 


61 


0.3 


38.5 


PG 0844+249 


158 


38 + 26 




2198 ± 174 


-1061 ± 175 


-558 + 51 








PG 1411+442 


-288 


-83 ± 18 




2362 ± 91 


-1265 + 92 


-709 ± 29 








PG 1448+273 


-106 


-38 + 31 


1288 


2129 ± 215 


-1102 ± 217 


-653 ± 67 








PG 1617+175 (U) -95 


-7 + 33 




2737 ± 234 


-1376 + 236 


-748 ± 69 










Median 




666 


2198 ± 174 


-1102 ± 175 


-653 ± 51 


64 


0.3 


38.5 



Note. — f The sum of the mass outflow rates from the multiple Gaussian components in the galaxy. Col.(l): Galaxy name. Col. (2): Heliocentric 
velocity of the absorption component as tabulated in Tables [3] and [5] Col. (3): With the exception of PG 1211 + 143 and PG 1440+356, this is the outflow 
velocity measured relative to the HI systemic velocity measured from HI emission using values tabulated in column (2) of Table [2] or [4"1 Vq, hi — Vq if no 
emission is detected. For PG 1211 + 143 and PG 1440+356, this velocity is relative to the systemic velocity measured from optical emission lines, since the 
detection of HI emission is uncertain (Figures [5] and [6j. Col. (4): Integrated HI column density based on the absorption model from Tables [3] and [5] Col. 
(5): Full-width at zero intensity based on the best-fit Gaussian absorption model and its associated error. Col. (6): The maximum outflow velocity where 
V raax = Vq, hi - ^FWZI. Col. (7): The velocity at which 98% of the gas has lower velocity where V 98 % = Vq, hi - 2( F ^™ ). Col. (8): Derived mass 
outflow rate calculated using Equation 5 an d assuming the o utflow is in a wide-angled cone where the mass is flowing into a solid angle of tt from a radius 
of 1 kpc following the same assumptions as Morganti ct al. (2005b). Col. (9): Same as column (8) but for a jet whose opening angle is assumed to be 5° 
launched at a distance of 0.5 kpc. Col. (10): Log of the kinetic energy outflow rate as calculated using Equation 6 for the jet scenario. 
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Fig. 1. — Full intensity GBT spectra of 25 galaxies whose velocities did not overlap with RFI signals. Columns from left to right are 
ordered by source type: ULIRGs — FIR-strong quasars - FIR-weak quasars. Each colu mn i s organized with increasing infrared luminosity 
going down. The targets with marginal or uncertain detections as noted in Section 13.21 are labeled with M and U, respectively, after 
their names. The x-axis is displayed as velocity relative to the heliocentric radial velocity of the galaxy given in Table [T] and identified 
by the vertical dashed line. The regions with hash marks are frequencies where the data are affected by persistent RFI spikes in some 
observing sessions, and the black regions are those where the data are affected by RFI spikes in all observing sessions. The best Gaussian 
fits to the absorption features are highlighted in blue, where the dotted lines are individual Gaussian components (when there are multiple 
overlapping components) and the solid curves are the sums of these components. The kinematics of these systems change from narrow 
absorption features likely due to outflows supported by the starbursts in ULIRGs, to broad absorption features from jet-driven outflows 
supported by the AGN in FIR-strong quasars, and finally to the emission profiles more characteristic of gas in rotation in the FIR-weak 
quasars. 
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Fig. 2. — The distribution of equivalent widths of the emission and absorption components as listed in Tables [2] and [3] The absorption 
components are shown as having negative EWs. Targets with only upper limits to their continuum flux densities are excluded from this 
plot. Along the merger sequence, the ULIRGs have the weakest features in both emission and absorption, the FIR-strong quasars have the 
strongest absorption features with moderate strength emission features, and the FIR-weak quasars are dominated by the strongest emission 
features. 
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Fig. 3. — The velocity widths of the components measured in emission (top; W50 in Table[2j and absorption (bottom; FWHM in Table[3]l. 
The colors represent ULIRGs (black), FIR-strong quasars (red), and FIR-weak quasars (blue). The width of the HI emission in the top panel 
may be inversely correlated with the SFR. The galaxies with the weakest star formation (FlR-weak quasars) have the strongest emission 
features, while the galaxies with the strongest star formation (ULIRGs) have the narrowest velocity widths. The velocity widths of the 
FIR-strong quasars fall between those of the ULIRGs and the FIR-weak quasars. There also seems to be a weak correlation between the 
FWHM in absorption and the SFR for ULIRGs, but this result is very uncertain as ULIRGs tend to have multiple absorption components 
for a single SFR. There is no apparent correlation between the SFR and the absorption parameters in the quasars, suggesting that the 
outflow-related absorption components are not the result of star formation. 
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Fig. 4. — The spectra of the 13 detected galaxies that show polarized HI spectra. The organization of this figure is the same as that of 
Figure[T] The colors represent the two different polarizations (red: XX; blue: YY). Polarized spectra are clearly seen in PG 0007+106, PG 
0050+124, PG 1119+120, PG 1211+143, PG 1351+640, PG 1411+442, PG 1440+356, and F15250+3608. Weak polarized features may 
be present in PG 0844+349, PG 1448+273, PG 1501+106, PG 1617+175 (U), and PG 2130+099. The FIR-strong quasars show the most 
distinct polarized features, perhaps reflecting the presence of polarized radio jets in these objects are contributing polarized continuum 
emission against which the absorption features are seen. 
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Fig. 5. — Left: GBT spectra of PG 1211+143 obtained at four epochs separated by at least one month, with the YY spectrum above and 
the XX spectrum below. The shifts in the central velocity of the deep absorption component are denoted by the vertical dotted lines in 
their corresponding colors. The central velocity shifted by ~120 km s over the seven-month period covered by our observations. Right: 
The 2011 November spectrum broken down into four 30-minute periods. There is clear variability in the shape of the spectrum, suggesting 
a varying co ntinuum a gainst which the absorption is detected. If the variability is real, it is consistent with the jet-induced outflow picture 
suggested bylMorganti (2011]). 
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Fig. 6. — Same as Figure[5] but for PG 1440+356. Short timescale variations are seen in this source, similar to those in PG 1211+143. 
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Fig. 7. — Background-subtracted 0.5—10 keV X-ray light curves of PG 1211+143 (left panel) over 6.5 years through four separate 
observations, and of PG 1440+356 (right panel) over one week through three separate pointings, from XMM-Newton. For each epoch, each 
point represents the EPIC-pn count rate in 1800 s intervals starting with the beginning of the observation in the first bin. The dotted lines 
are the median rates for each specific epoch. The periods of high background flares have been removed from the plot and the calculations 
of the median rate. The error bars are 3— a. The AGN have both high and low flux states and are clearly variable on the 30-minute time 
scale seen in Figures [5] and [6] 
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Fig. 8. — The HI screen and source configuration that explains the variable polarization-dependent spectra observed in PG 1211+143 
and PG 1440+356. In this picture, the HI screen is extended, covering both the core and the radio jets. The jets produce the polarized 
continuum against which the absorption features are detected. As the central source fluctuates, the ratio of central-to-jet emission also 
changes, resulting in the observed spectra. 
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Fig. 9. — Complex decomposition of the polarized HI spectra assuming that the difference in the two polarizations is due to more HI 
absorption in one polarization than the other. For each galaxy, the upper panel is the difference spectrum between the two polarizations, 
i.e., the spectrum in the polarization with the more negative flux density (Poll) minus the spectrum in the polarization with the higher 
flux density (P0I2), plotted in purple. The black curve shows the Gaussian model component(s) best fitting to the difference spectrum as 
tabulated in Table [5] The lower panel shows the residuals of the best fit model as applied to the Poll spectrum and the P0I2 spectrum to 
demonstrate the quality of the fit. The colors in all of the bottom panels, the x-axis, and the hashed rectangles are the same as those in 
Figure [4] The polarization dependence of the HI spectra can be explained by broad HI absorption affecting one polarization more than the 
other. 
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Fig. 10. — Same as Figure [2] apart from the use of the complex fitting results when available (Tables l4l and l5t . The general trends with 
stage along the merger sequence are in agreement with those in Figure [2] 
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Fig. 11. — Same as Figure \3\ apart from the use of the complex fitting results in Section 15.21 when available (Tables [4] and [5j . The general 
trends are in agreement with those in Figure [3] 
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Fig. 12. — The distribution of FWZI of the absorption components in galaxies with detected outflows with Vq hi l ess than -50 km s 1 
(see Table There is a tentative progression of velocity widths from ULIRGs to FIR-strong quasars to FIR-weak quasars. 



